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Abstract 

A symbolic computational algorithm which detects " linear " ' solutions of non- 
linear polynomial differential equations of single functions, is developed through- 
out this paper. 

1 Introduction 

The problem of obtaining general solutions of differential equations via symbolic algo- 
rithms has been studied in the past by many authors, [21)0],|Hlj|ni)[Zlj|H|- These algo- 
rithms allowed new calculation techniques to be accomplished, much more efficiently, 
faster and without approximation errors. 
In this paper we treat with differential equations of the form: 

p(x,y(x),y'(x),y"(x),...,y( n \x)) = (1) 

where p is a polynomial function and y{x) a real function of a single variable. Our 
aim is to discover possible " linear " solutions of By the term "linear" we mean 
solutions which can be obtained by solving linear differential equations. Our approach 
is focused on the construction of an algorithm which faces the problem symbolically. 
What this algorithm is essentially doing is that helps us to rewrite p as follows: 

P = a(x, y',..., y^, WV^)[W , fe + W hk y + W 2 , k y' + ■■■ + y {k) ] joA 

[[Wo, fc+M + Wi.fc+i.11/ + W 2 , k +i,iy' + ■■■ + y {k) ]'] hA 
■ ■ ■ [[Wb,„,i + W ltn>1 y + W 2 , n ,iy' + ■■■ + j (t) ] ( "~ tl ] ) "- w 
+ --- + c lJ (x,y',...,y^- 1 \W t ^ v ) 
[W , k + W ltk y + W 2 , k y' + ■■■ + y^Y"- 1 
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[[Wo.k+i.u + wi.k+1,*!/ + w 2 ,k+i,,y' + ■■■ + ?/ fe) ]T' 1 -" 

• ■ ■ [\W ,n,v + W hn , v y + W 2 , n , u y' + ■■■ + y ( fc )](»- fc )]i»-*.- + i? 

Where k is given, n is the order of p, and j a ^ are specific positive whole numbers. The 
quantities Wi t(TtV are undetermined parameters which can take certain values, Cj are the 
coefficients, depending from the parameters Wi^ a>v and R a polynomial of the variables 
x, y,y', . . . , and the parameters Wj,o-, y called the remainder. Afterwards, we seek 

for those values of the parameters which eliminate the remainder. If this is possible, 
then the linear differential polynomial Wo,k + Wi t kU+ • ■ ■ +jr fc_1 \ with Wb,fcj Wi,fe) ■ ■ ■ 
evaluated over those values which vanish the remainder, is a factor of p, (where the 
operation of differentiation has been taken under consideration) . This means that any 
solution of Wo t k + Wi^y + • ■ • + j/'* -1 ' = is a solution of the equation QJ, too. Since 
we do it for every k = 0, . . . ,n, we collect likewise, all the " linear " solutions of the 
equation p = 0. 

Our method is an extension of a similar procedure, introduced by the author and 
applied in the study of difference equations and feedback design Its main merit is 

its computational orientation. It turns to be a useful tool, implemented on a computer 
machine and gives useful results. Moreover, despite our method resembles with the 
approach of differential algebra, 0, there are some differences. Concretely, (i) We 
are working with a single differential polynomial whilst Ritt's algorithm deals with 
sets of differential polynomials, (ii) The existence of the parameters Wi^. v permits us 
to find classes of linear solutions. We can then select among them, these particular 
solutions which satisfy additional conditions. For instance, we can search for those 
values of Wi i(Ti¥ ,, if any, which do not only eliminate the remainder but also yield stable 
linear factors, (iii) In the classical Ritt's approach we find the minimum number of 
differential polynomials which generate a differential ideal. In our method we check 
if a given polynomial belongs to a differential ideal, produced by linear differential 
polynomials. Throughout the text, R and Z + will denote the sets of real numbers and 
positive integers, correspondingly. 

2 The Algebraic Framework 

Let R[x] be the ring of polynomials of a single variable with real coefficients. This 
polynomial ring is a differential ring too, with the usual derivation 0,0 • Let y(x) be 
a real function and y^(x), i = 0, 1, 2, . . . its derivatives. A differential polynomial p, in 
y(x) or shortly in y, is a polynomial in y and its derivatives with coefficients in H[x]. 
p can be written as follows: 

if n 

where s\ € R and some of the exponents a\,9i t \ € Z + are not equal to zero. The 
number n, which represents the highest order derivative of y(x), is called the order 
of p. An equation of the form p = 0, with y(x) as unknown function, is called a 
polynomial differential equation. Any function which satisfies it, is called a solution or 
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a general solution. An expression of the form: L = Yl7=o a iy^{ x )i a i € R, is called a 
linear differential polynomial and the equation I = 0a linear differential equation. Its 
solutions are called " linear " solutions of order n. 

The highest derivative of y(x), appeared in the polynomial p, is called the leader. Let 
Pi = six ai nr=o[y^H a; )] ei ' 1 7 -P2 = s 2 x a2 Yli=o[y® i x )] 0i ' 2 be two, not identical, terms of 
p. This means that there is at least one index k, 1 < k < n, such that 9k,i ^ Ok, 2 or 
= 9i,2> i = 0, 1, . . . n and a\ ^ a 2 . We say that the term p 2 is ordered higher than 
Pi with respect to lexicographical order and we write p\ -< P2, if either there is an index 
s such that # Sj i < 9 S ,2 and 9j\ = 9j t 2, j = s + 1, . . . , n or 6^1 = 6^2, t = 0, 1, . . . , n and 
a± < <22 • By means of this rank we can order all the terms of p in an ascending way. 
The term which is ordered higher, is called the maximimum term of p. 
The differential ideal, generated by a finite set of differential polynomials: 
<I> = {$1, $2, • ■ ■ , ^m} and denoted by [$] is a set which consists of all differential 
polynomials that can be formed of elements in 3? by multiplication with arbitrary 
polynomials, addition and differentiation. 

Let W = {Wi t a yV } be a set of undetermined parameters, taking values in R. A Formal- 
k- Factorization of p, denoted by Formal(p, k), is an expression of p of the form: 

V 

Formal{p, k) = } j c u (x, y'(x), . . . , y^^ix), W iia , v )- 

n — k 

■ II [(Wb.fc+i,,! + W hk+t ^y(x) + W 2tk+itli y'(x) + ■■■ + y {k) + R (2) 

i=0 

where the coefficients c M (a;, y'(x), . . . ,y^ k ~ r) (x), Wi^ tV ) and the remainder R are poly- 
nomials of the terms x, y'(x), . . . , y( k ~^ (x), Wi ta - ttp only. Some of the exponents ji € Z + 
may be equal to zero. Sometimes, © is written briefly as 

v n — k 

Formal(p, k) = ^ c M J[ + i? 

where we used the notation L r ^ for the linear differential polynomial Wo^+i,^ + 
Wi,k+i, M y(x) + W 2 ,k+-^y'{x) +••'• + y( fe ' (a). 

We can take different expressions of the Formal(p, k) of a concrete differential poly- 
nomial p, by giving to the parameters Wj ](7i¥ , certain values. Such procedures are 
called evaluations of the Formal (p, fc). A most rigorous approach is the following: Let 
W = {Wi t<T>{ p} be the set of the variables, appeared in the Formal - k - Factorization of 
a given polynomial p. By arranging the parameters in an increasing order we form the 
vector W = (Wi h ^ h , Vh )h=i,2,...,7i- Let r = (ah)h=i,2,...,n be a vector of real numbers, 
which has the same length with the vector W. We say that the parameters W follow 
the rules r and we write W — *• r if the following substitution is valid: Wi h _ ah _ iph = ah, 
h = 1, 2, ... ,ri. Let M a set of rules, M = {ri, ■ ■ ■ , ta} then 

A v n — k 

Formal(p, k) \ M = (J c M . . [] [(W„ 

u=l u=l i=0 
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+W 2 M +l ^y'{x) + --- + y^(x)f^ +R with W - r„ £ M} 

The set of substitutions M, may be finite or infinite. When we evaluate Formal(p, k) 
over M, the linear differential polynomials Lj,^ and the remainder i? take specific 
values, these are denoted by | and i? | A case of particular interest is when 

we can find values of the parameters which eliminate the remainder R. Whenever this 
happens, p is a " combination " of linear differential polynomials or, in a more formal 
language, p is a member of the differential ideal produced by these linear differential 
polynomials. Relevant is the following theorem: 

THEOREM 1. Let p be a differential polynomial of order n, let k be given and 

Formal(p, k) = X^=i c m 0^=0° [-^i*^*]" 7 "* + ^ its Formal-k- Factorization. Let us suppose 
that there is a set of rules, denoted by TZ, which eliminates the remainder R, i.e. 
R \ n = 0, then p e [ \ n ,i = 0, . . . , n - k, /j, = 1, . . . , v\. 

PROOF. The proof comes straightforward from the definition of the Formal-k- 
Factorization. □ 

It is obvious that if the linear differential equations Lj ;/J | = have a common 
solution, this is a solution of the nonlinear equation p = 0, too. This is the cornerstone 
of our approach. 

EXAMPLE 1. Let us consider the differential polynomial 

p = Av"'-A( v "f + y'v"-^ > {y'?-i 

then, 

Formal (p,2) = 4(Wo,3,i + Wi,3,iy+W2,3,iy / +y / / - 4 (^o,2,2 + Wi,2,2y+W 2 ,2,2y / +y / ') 2 

+ [8W ,2,2 - 4W 2 ,3,1 + 8VFi, 2 ,2 + (8VF 2 ,2,2 + 1)?/] • (VF ,2,3 + Wl,2,32/ + ^2,2,32/' + 2/") + # 

where the remainder R is R — 4Wo 2 2 _ 8Wo, 2,3^0,2, 2 

+ 4^0,2,3^2,3,1 ~! +(8T^l,2,2Wo,2,2'-8VFi,2,3VFo,2,2-8VF ,2,3^1,2,2+4T4 / i,2,3W / 2,3,l)y 

+ ( 4W/ l 2 2.2-8VFi,2,2VFi, 2 ,3)y 2 +(8^2,2,2^,2,2-8^2,2,3^0,2,2 -^,2,3-4^1,3,1-8^,2,3^2,2,2 

+4^3^2,3,1)^ +(4^ 

— 8Wi,2,3W2,2,2 —8^1,2,2^2,2,3)2/2/'- The following rules eliminate the remainder: ri = 
{Wo, 3,1 = s, VFo,2,2 = w, W2,2,3 = ¥?, VF ,2,3 = k, VFi, 3 ,i = \u- (4 " 4 fc 1)v:> , Wi, 2 ,3 = 0, 
W 2 , 2 .2 = -§, Wi.2,2 = 0, W2.3.1 - 2w + and r 2 = {VF , 3 .i = a, W . 2 ,2 = u, 

W 2 ,2, 3 = <p, VFo, 2 ,3 = k, VFi, 3 ,i = ±[p + (w - 2fc)(fc - 4(w - 2fc)p)], VFi, 2 , 3 = 0, 

VF 2 ,2, 2 = 2ip + i, Wi )2 ,2 = 0, VF 2 , 3 ,i = 2u + ^|^} with w,<p,k,s E R. Indeed, for 
instance 

Formalip, 2) | rj = 4( S - 1 ( w + 2/ + + 2/' + y")' 

-4( W - | + 2/") 2 + - ^) • (" + + 2/") 

and the differential ideal which contains is [s — \ (uj + ( 4 " y, + + x ~%£ ^ y' + 

y" , uj — ^ + y" , uj + ipy' +y"] u, <p, s, k, € R. By setting w = i^ = s = 0,fe = lwe take the 
simplified ideal [y", —\ + y"] ■ We can take analogous expression for Formal(p,2) | . 
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3 Detection of the Linear Solutions 



The scope of this section is to present the algorithm which constructs for a given k, a 
special Formal(p, k) with a linear differential polynomial as a common factor to each 
term. We denote this factor by L c ,fc- Afterwards, by finding proper sets of values for 
the parameters Wj jCr ,y, we eliminate the remainder. It is then clear that any solution 
of the linear equation L c .k = 0, where the polynomial L c .k has been evaluated over 
this set, is a solution of the original system, too. By repeating the whole procedure 
for every k = 0, . . . ,n we discover all the linear solutions. As we pointed out, the 
crucial issue is how can we eliminate the remainder. This is carried out by solving a 
system of algebraic equations. Finally, we have to elucidate that in this paper we do 
not take into account initial conditions. We are only focused on how we obtain general 
solutions, that is solutions which " contain " constants. We present now the algorithm 
upon discussion. 

Let us suppose that an algorithm which solves an algebraic system of polynomial equa- 
tions, is available. These algorithms are classical in computational algebra and there 
are many of them in the literature joj- We name a such algorithm as SysAlgEqs. 

THE DIF-FORMAL ALGORITHM 

Input: 

• A differential polynomial p of order n. 

• A set of undetermined parameters W = {Wi, a ,ip\, taking values in R. 

Output: The quantities S k , k = 0, . . . , n 
FOR k = TO n 
Step 1: R — p, fi — 0. 

Step 2: REPEAT the following steps UNTIL R does not contain terms of order > k. 
Step 2a: Set fj, = fi + 1, 

Step 2b: Find the maximum term of R, with respect to the lexicographical order. 
We denote it by 

n 

pp = 8 II -x*> -n[y (i) (*)]*'■" 

where a M ,Ai i(J are positive integers and y^\x) the derivatives of y(x) of order 
i. At the first iteration s M is a constant, then it becomes a function of the free 
parameters Wi t<y , v as well. 

Step 2c: Construct the linear formal differential polynomials: 

L c , k = W , h + Wi, k y(x) + W 2 ,kV'(x) + ■■■ + W k , k y (k - 1} (x) + y (k) (x) = 

k-l 

= W , k + J2w j+ i, k y U \x) + y ( - k \x) 

3 = 

k-l 

Li, k = W ,k+i,» +^2Wj + i,k+i,»y U) (x) + y (k] '(a;), i = l,...,n-k 

3=0 
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Step 2d: Execute the operation: 

k— 1 n—k 
i=0 i=l 

END of REPEAT 

Step 3: By means of the SysAlgEqs-Algorithm we find the set Sk of those values of the 
parameters, Wi , CT ,w, which eliminate the remainder. In other words: R I =0. 

' ' 1 ^k 

END of FOR 

It is obvious that the DIF-FORMAL algorithm terminates after a finite number of 
iterations. 

THEOREM 2. Let S k , k = 0, . . . , n be the outputs of the DIF-FORMAL Algorithm. 
If Sk for some values of fc, then the solutions of the linear differential equations 
L c k L =0 arc solutions of the nonlinear polynomial differential equation p = 0, 
too. 

PROOF. Let p a differential polynomial and k fixed. By substituting backwards 
the successive results of the step 2d we find that 

m fc— 1 n—k 

p = e v=°" n • • n a - 

^i — 1 Z— 2 — 1 

This is the Formal-k- Factorization of p with c M = s^x "* YiiZoly^]^''* an d ^c,fc as a 
common factor in all the terms but the remainder. We denote it by CFormal(p, fc). Let 
us now suppose that Sk ^ 0, this means that R L =0 and thus, the linear differen- 
tial polynomial £ c evaluated over S^, is a factor of every term of CFormal(p,k) I = 

1 <^k 

0. This implies that any solution of the linear differential equation L c k = is 

' ' Sk 

also a solution of the nonlinear equation p = 0. Since this argument is true for any fc, 
the theorem has been proved. □ 

The above result can be restated, using ideals, in the following way. 

COROLLARY 1. Let S k , fc = 0, . . . , n be the outputs of the DIF-FORMAL Algo- 
rithm. If Sk ^ for some values of fc, then p £ [ L c ,k L ], fc = 0, . . . , n. 

' ^k 

In order to eliminate the remainder we have to solve a system of algebraic polyno- 
mial equations. This can be done via several methods. Groebner basis, 3 , is a popular 
powerfull tool, with satisfactory results. 

We can extend the algorithm posted above toward different directions, (a) Instead 
of the constant term in the common linear differential polynomial we can use a single 
polynomial of a: of a given degree with parametrical coefficients. These coefficients can 
be determined through the remainder elimination, too. (b) Instead of linear differen- 
tial polynomial, with constant coefficients, we can use linear differential polynomials 
with variable coefficients which correspond to solvable differential equations (Euler 
equations, for instance). The last two cases are currently under studying. 
EXAMPLE 2. Let us consider the differential equation 

{y"f - 2y'(y") 2 - %") 2 + {y'fy'" + W" + V" = 
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or p = 0, where y = y(x) is the unknown function and the order of p is n = 3. 
We want to find all the "linear" solutions which are included into this equation. The 
application of the DIF-FORMAL algorithm gave the following results: y(x) =c,ceR, 
y(x) = Ci + xc 2 , y(x) = c\e x — 2x + c 2 , Ci, c 2 G R. To clarify our ideas and to indicate 
how the algorithm works in practice, we shall present the case k = 2 in details. Since 
k = 2 we are going to detect linear polynomials of second order included into the 
original equation. At the first iteration the maximum term is [{y') 2 ]y"' ■ We construct 
the linear polynomial: 

L c ,2 = Wo, 2 + W h2 y + W 2 , 2 y' + y" 

and we execute the subtraction: pi = p — [{y') 2 ] ■ L' c2 . This operation will eliminate 
the y'" term. In the next iteration Ay'y 1 " is the maximum term. Since we are looking 
for a common second order linear equation included into p, we use the same L c2 and 
we calculate p 2 = P\ — %' • L Cj2 . Working this way we finally get: 

Formal(p, 2) = [{y' f + Ay' + 4] • L' c 2 + L% 2 + [-3W , 2 - SW h2 y- 

~3W 2 , 2 y' - 2y' - A\L 2 cfi + [iW 2 a + 6W h2 W , 2 y + 6W 2 , 2 W , 2 y' + AW Q , 2 y'+ 

+8VFo, 2 + 3W 2 2 y 2 + SW 2 ^ 1 ) 2 + 3W 2 , 2 (y') 2 - 4W 2 , 2 + 8W h2 y + 4W 2 , 2 y'+ 

+ / LW h2 yy' + Wi >2 W 2>2 yy']L ca + R 

where R = (-3^i. 2 W 2 2 -Wi, 2 Wo. 2 +4W h2 W 2t2 )y +(-3W , 2 W? 2 -4W? 2 )y 2 -W? 2 y 3 
+(-3W 2 ,2Wo 2 2 -2Wg 2 -4W 2 ' 2 W , 2 +4Wi' 2 -m x . 2 )y' +(-4W ' 2 W lfi 
-6W 2 W 2 2 W X 2 -4W 2 2 Wi 2 )yy' +(-3W 22 W? 2 -2Wl 2 )y 2 y' +{~3W 2 W% 2 -3W 2 W 2 2 
-4VF li2 )(y') 2 +(-3Wi,2^ 2 2 2 -3W h2 W 2 , 2 )y(y r ) 2 +(-Wl 2 -W 2 2 -W h2 )(y r f +AW 2j2 W^ 2 
-W 3 >2 - 4W 2 2 . The values {W oa = -2, Wi, 2 = 0, W 2a = -1) and {W oa = 0, Wi, 2 = 
0, W 2 ,2 = 0), eliminate the remainder. By substituting them to L c . 2 we get the linear 
differential equations y" — y' — 2 = 0, y" = 0. The general solutions of those equa- 
tions will be solutions of the original nonlinear differential equation too. Therefore 
y(x) = c\e x + c 2 — 2x, y(x) = c 2 x + c\ are solutions of the equation p = 0, too. Work- 
ing similarly, we obtain for k = and k = 1 analogous results. What we have actually 
proved is that p is a member of the differential ideal [y" — y — 2, y"]. 
EXAMPLE 3. Let us illustrate now the case where, instead of constant terms, we have 

polynomials of x. We consider the differential equation: 

{x 2 - x)y' + xy" - x 2 y"' + (-2x 3 - 3x 2 + 3x) = 

In this case we work with the extended algorithm which uses common linear factors of 
the form: y' + W 1A y + W- 2fi .xx 2 + W- lfi ,iX + W ,o,i- This modified DIF-FORMAL 
Algorithm will give that Wi t i — — 1, W- 2 .o,i = 1, W-1,0,1 = 3 and Wo,o,i free or 
Wi.i = 0, W- 2 . o,i = 0, VF_i,o,i = —2 and Wo,o,i = —5. The corresponding general 
linear solution is y(x) = c\e x + x 2 + 5x + c 2 . 
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